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I. INTRODUCTION 

Observations of charmless three-body baryonic B de- 
cays have been reported recently by both the Belle and 
BaBAR collaborations • A common feature of these 

decay modes is the peaking of the baryon-antibaryon 
mass spectrum near threshold. This feature has stim- 
ulated considerable interest among theorists as a key ele- 
ment in the explanation of the unexpectedly high branch- 
ing fractions for these decays 0, [H] . 

In the standard model, the B^ Apu^ decay proceeds 
through tree level & — > m and penguin b ^ s amplitudes. 
It is of interest to study the structure of the decay ampli- 
tude in the Dalitz plane to test theoretical expectations. 
The weak decay A — > pTr"*" is spin self-analyzing. Since 
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the s quark carries the A spin, the V-A transition & — > s 
leads to the expectation that the A is fully longitudinally 
right-polarized at large A energy in the -B" rest frame [a] ■ 
This channel may also be used to search for direct CP 
violation. 



II. DATASET AND SELECTION 

The data sample consists of 467 x 10^ BB pairs, corre- 
sponding to an integrated luminosity of 426 fb^^, col- 
lected at the T(4S') resonance with the BaBAR detec- 
tor. The detector is described in detail elsewhere [7|. 
Charged-particle trajectories are measured in a tracking 
system consisting of a five-layer double-sided silicon ver- 
tex tracker (SVT) and a 40-layer central drift chamber 
(DCH), both operating in a 1.5-T axial magnetic field. 
A ring-imaging Cherenkov detector (DIRC) is used for 
charged-particle identification. A CsI(Tl) electromag- 
netic calorimeter (EMC) is used to detect and identify 
photons and electrons, while muons and hadrons are 
identified in the instrumented flux return of the magnet 
(IFR). A BaBAR detector Monte Carlo simulation based 
on GEANT4 [8] is used to optimize selection criteria and 
determine selection efficiencies. 

We reconstruct A candidates in the A ^ pir decay 
mode as combinations of oppositely charged tracks, as- 
sign the proton and pion mass hypotheses, and fit to a 
common vertex Q. Combinations with invariant mass 
in the range 1.111 - 1.121 GeV/c^ are refit requiring the 
track pairs to originate from a common vertex and con- 
straining the mass to the world-average A mass po| . 
Candidate mesons are formed by combining A candi- 
dates with two additional oppositely charged tracks, each 
with momentum transverse to the beam greater than 
50 MeV/c. 

Measurements of the average energy loss (dE/dx) in 
the tracking devices, the angle of the Cherenkov cone 
in the DIRC, and energies deposited in the EMC and 
IFR are combined to give a likelihood estimator La for 
a track to be consistent with a given particle hypothe- 
sis a. We require that the yl-decay proton candidates 
satisfy the particle-identification criteria Lp/Lx > 0.33 
and Lp/L^^ > 1 to discriminate from kaons and pions, 
respectively. The candidate protons, which are assumed 
to originate from the B^ decay vertex, are analyzed with 
a selection algorithm based on bagged decision trees [ll| 
which provide efficient particle discrimination, retaining 
96.4% of the signal candidates and 17.8% of the back- 
ground. The candidate pions from the B^ vertex are 
required to pass a similar selection algorithm, tuned to 
discriminate pions, that retains 98.8% of the signal and 
66.8% of the background. A Kalman fit [1^] to the full 
decay sequence is used to reconstruct the B° vertex us- 
ing the position of the beam spot and the total beam 
energy as kinematic constraints. Only candidates with 
a fit probability Pvtx > 10~^ are considered, a require- 
ment that retains 94.4% of the signal and 58.2% of the 



background. 

The primary background arises from light-quark con- 
tinuum events e"*'e~ ^ il — u,d,s,c), which are 
characterized by coUimation of final-state particles with 
respect to the quark direction, in contrast to the more 
spherical BB events. Exploiting this shape difference, 
we increase the signal significance using event-shape vari- 
ables computed from the center-of-mass (CM) momenta 
of charged and neutral particles in the event. For each 
event, we combine the sphericity [l^l, the angle between 
the B° thrust axis and detector longitudinal axis, and 
the zeroth and second-order Legendre polynomial mo- 
ments [l3] of the tracks not associated with the recon- 
structed B candidate, into a Fisher discriminant p^ . 
where the coefficients are chosen to optimize the sepa- 
ration between signal and continuum-background Monte 
Carlo samples. We find that the selection using the op- 
timal cut on the Fisher discriminant retains 72% of the 
candidates from the signal Monte Carlo sample and 8% 
from the continuum-background Monte Carlo sample. 

To further reduce the combinatoric background, we 
take advantage of the long mean lifetime of A particles 
and require that the separation of the A and B'~' vertices, 
divided by its measurement error, computed on a per- 
candidate basis by the fit procedure, exceeds 20. This 
criterion is optimized on Monte Carlo events and is ef- 
fective in rejecting 42% of combinatoric background that 
survives all other selection requirements, while retaining 
90% of the signal candidates. The only sizable B" back- 
ground is from the process B'^ Ac p Apir" , which 
we suppress by removing candidates with an invariant 
mass m{ATT~) within 5 standard deviations (20MeV/c^) 
of the nominal Ac mass [lO]. 

The kinematic constraints on i?" mesons produced 
at the T{4:S) allow further background discrimination 
from the variables toes a-nd AE. We define toes — 

\J (f +Pi -PbY lEf - ps^, where {Ei,pi) is the four 
momentum of the initial e'^e~ system and pb is the mo- 
mentum of the reconstructed B^ candidate, both mea- 
sured in the laboratory frame, and s is the square of the 
total energy in the e+e" center-of-mass frame. We define 
AE = E*^-^, where E% is the B° energy in the e+e" 
center-of-mass frame. Signal candidates have toes close 
to the B'^ mass and AE near zero. Candidates satisfying 
|A£:| < 100 MeV and 5.20 < toes < 5.29GeV/c2 are used 
in the maximum-likelihood fitting process. 



III. BRANCHING FRACTION 

We measure the branching fraction with a maximum- 
likelihood fit on the TOEs-Ai? observables of recon- 
structed B° candidates. The s'Plot technique [lB| is then 
used to determine the m{Ap) distribution and, after cor- 
recting for the nonuniform reconstruction efficiency, mea- 
sure the TO(ylp)-dependent differential branching fraction. 

We consider as signal candidates only reconstructed 5° 
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candidates in which all particles are correctly assigned in 
the decay chain. By self-cross-feed, we refer to events in 
which mesons decay to Apn and are reconstructed as 
signal candidates in which one or more particles are not 
correctly assigned in the decay chain. An example of such 
a misreconstruction is where the protons from the signal 
and A decays are interchanged. We define the prob- 
ability density function (PDF) in the A£'-toes plane as 
the sum of signal, self-cross-feed, and background com- 
ponents. The likelihood function is given by 



1 



N 



Y[{NsPs{ye) + NscfPscfiye) 
e=l 

+NBVB{ye)}, (1) 



where j/e = (wESei^^e), the product is over the N 
fitted candidates with Ns and Nb representing the num- 
bers of signal and background events, and Nscf = Nsfsd 
representing the self-cross-feed contribution. The three 
V functions are taken as products of one-dimensional AE 
and rngg PDFs. We are justified in this simplification by 
the small correlation between these two variables in our 
Monte Carlo sample. The toes PDF is taken as a sum of 
two Gaussians for the signal and an ARGUS function 
for the background. The AE PDF is taken as a sum of 
two Gaussians for the signal and a first-order polynomial 
for the background. Finally, the self-cross-feed contribu- 
tion shows a peaking component that is modeled as the 
product of a sum of two Gaussians in AE, and a single 
Gaussian in m^s. We determine fsd = 0.006 and the 
other parameters that characterize this background from 
fits to simulated events. 

We fit the means of the narrow AE and m^s signal 
Gaussians, the coefficient in the exponential of the Argus 
function, the linear coefficient of the AE background dis- 
tribution, and the event yields Ns and Nb- The means 
of the wide Gaussians are determined by applying Monte 
Carlo-determined offsets to the means of the narrow ones, 
such that only an overall shift of the fixed PDF shape is 
allowed. All other parameters used in the likelihood def- 
inition are fixed to values determined from fits to Monte 
Carlo-simulated events. 

Once the maximum-likelihood fit provides the best es- 
timates of the PDF parameters, we use the ^^lot tech- 
nique to reconstruct the efficiency-corrected m{Ap) dis- 
tribution and measure the branching fraction. The PDF 
is used to compute the s-weight for the nth component 
of event e as 



sVniVe) 



(2) 



where the indices n,j, and k run over the ric = 3 signal, 
background, and self-cross-feed components. The sym- 
bol V„j is the covariance matrix of the event yields as 
measured from the fit to the data sample. An important 
property of the ^'Plot is that the sum of the s-weights 
for the signal or background component equals the cor- 
responding number of fitted signal or background events. 



We have demonstrated with simulated experiments that 
the s'Plot is an unbiased and nearly optimal estimator 
of the m(Ap) distribution. To retrieve the efficiency- 
corrected number of signal events in a given m{Ap) bin 
J we use the s-weight sum 



eGJ 



(3) 



where the per-event efficiency e(xe) depends on the posi- 
tion = (rrij^p, cos 9 in the square Dalitz plane. Here 
9j[^ is the angle between the momenta of the pion and 
the A candidate in the Ap rest frame, and the efficiency is 
determined over a 20 x 20 grid in the Dalitz plane, using 
fully reconstructed signal-Monte Carlo events. The error 
o'[Ns,j] in Nsj is given by 



(4) 



An estimate of the efficiency-corrected number of signal 
events in the sample is given by the sum of the efficiency- 
corrected s-weights, or 



and the branching fraction is obtained from 



B {B° ^ Apn-) = 



Ns 



N 



BB 



B{A^ pir) ' 



(5) 



(6) 



where Nbb is the total number of BB pairs and 
B{A-^Ptt) = 0.639 ± 0.005 [lO]. Using a collection of 
Monte Carlo pseudoexperiments, in which signal candi- 
dates, generated and reconstructed with a complete de- 
tector simulation, were mixed with background candi- 
dates, generated according to the background PDF, we 
confirm that this procedure provides a measurement of 
the branching fraction with negligible biases and accurate 
errors. 

We can measure the CP-violating branching-fraction 
asymmetry by tagging the flavor of the (B^) meson 
with the charge of its daughter proton (antiproton). We 
repeat the maximum-likelihood fit described above in- 
cluding the partial rate asymmetry 



A = 



B{BO ApTT+) - B{B° Apn-) 
B{B" ApTr+) + B{BO Apn-) 



(7) 



as a free parameter. We reduce the effect of systematic 
differences in particle-identification efficiencies between 
protons and antiprotons, and between positive and neg- 
ative pions, by performing the fit on a sample of re- 
constructed candidates, where protons and pions that 
originate from the A decay satisfy the same particle- 
identification criteria as those imposed on the protons 
and pions that originate from the vertex. 
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IV. A POLARIZATION MEASUREMENT 

We study the three orthogonal components of the po- 
larization of A candidates reconstructed in the 5° — » 
ApTT~ decay as a function of E*^, the A energy in the B° 
rest frame @. The distribution of the helicity angle 6'h 
for the A decay is given by 



1 dT _ 1 |- 
f dcos On ^ 2 ^ 



iP (S*-)cos0h] 



(8) 



where 0h is the angle between the antiproton direction, 
in the A rest frame, and either (1) L, the unit vector 
in the direction of the A in the rest frame; (2) T, 
the unit vector along the direction of the cross product 
between the momenta, in the B'^ rest frame, of the proton 
and the A; or (3) TV = L x f . The symbol P {E\) is 
the component of the A polarization in the L, T, or N 
direction as a function of and a^i is the A decay- 
asymmetry parameter [l^. CP conservation in 5° — > 
Ap'K~ decays implies that 

'^aP[LM],B°^APtt {E*a) = 01aP[L,N],B°->Apt! {E*a) ' (9) 

while the product olaPt changes sign under CP con- 
jugation. We use these relations to fit the B^ and 
candidate samples together. 

We use a maximum-likelihood fit in mEs, Ai?, E*^, and 
cos to measure the polarization as a function of E*^ 
along each of the three axes defined above. We divide the 
E*j^ range into three bins with boundaries 1.10, 1.53, 1.80, 
and 2.40 GeV, chosen in order to have similar numbers 
of signal events in each bin. We define a PDF as the sum 
of signal and background components. The likelihood is 



-{Nk.s+Nk.B) 



k=l 



- Nk,BV'k.B{Ze)VB{ye)], (10) 



where we have divided the observables into two sets 
Ue = (mEs, AS) and Ze = {cos6h, E'^), and the products 
are over the three bins in E*^ and over the N/. events 
that populate the kth bin, where Nk^s and Nk,B repre- 
sent the numbers of fitted signal and background events. 
The Vs^BiVe) PDFs are the same functions used in the 
branching- fraction measurement. However the self-cross- 
feed component is not included since it corresponds to a 
negligible fraction of the signal events. For the kth bin in 
E'j^, the signal (cos 6'h, -E^) PDF is written as the product 
of the differential branching fraction of Eq. O times the 
signal-reconstruction efficiency e (6'h , E*^ : 



V'k.s (0H,i?l) - [l + {c 



iPjfcCOSfe] (11) 



where the {a^P^k are fit parameters. The signal- 
selection efficiency is measured with a sample of recon- 
structed signal-Monte Carlo events that pass the same 



CD 

o 




1.2 1.4 1.6 1 



(GeV) 



FIG. 1: Reconstruction efficiency measured on the Monte 
Carlo signal sample, as a function of cos^h and E\. In this 
plot On is the angle between the antiproton direction, in the 
A rest frame, and L, the unit vector in the direction of the A 
in the rest frame. 



selection criteria as those used to define the data sam- 
ple. We bin the signal efficiency in 20 x 20 rectangu- 
lar boxes that cover the allowed region of the i?^-cos 6'h 
plane (Fig.[T|). 

The background 6'h distribution is modeled as a lin- 
ear combination of Chebyshev polynomials up to fourth 
order. The four coefficients that define the linear com- 
bination are fitted independently for each of the three 
bins in E*^. We study the 6-a distribution of back- 
ground events using candidates in the sideband region 
Toes < 5.27GeV/c^, and find it to be nearly independent 
of Toes- We consider this insensitivity as an indication 
that the shape of the background 6'h distribution is the 
same for events in and out of the signal region. 

We have confirmed that this PDF representation does 
not bias the polarization measurement by performing 
pseudoexperiments in which signal candidates, gener- 
ated and reconstructed with a complete detector simu- 
lation, were mixed with background candidates gener- 
ated according to the observed helicity distribution in 
the Toes < 5.27GeV/c^ sideband. The number of sig- 
nal and background candidates are chosen to match the 
characteristics of the data. 



V. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties in the branching-fraction 
measurement are listed in Table I and classified as overall 
uncertainties, uncertainties associated with event selec- 
tion, and uncertainties associated with fitting the event 
distribution. We study the uncertainty due to tracking 
efficiency by comparing data and Monte Carlo for a sam- 
ple of T-pair events, in which one r decays to one charged 
track and the other r decays to three charged tracks. 
We separately study the tracking efficiency of A decay 
products using an inclusive sample of yl — s- p7f candi- 
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TABLE I: Systematic uncertainties on the branching-fraction 
measurement. "Total" is the sum in quadrature of all the 
individual contributions. 





Source 


T T j_ ■ j_ /fry \ 

Uncertamty (%) 


Overall 


Tracking efficiency 


2.4 




PID efficiency 


1.4 




MC statistics 


0.4 




BB counting 


1.1 




B B"/BB traction 


3.2 




A ^ pn branching fraction 


0.8 


Event selection 


Event shape 


1.0 


requirements 


Fit probability 


1.0 




A fiight length 


2.8 




A mass 


2.4 




Ac veto 


0.5 


Fit procedure 


Likelihood parameters 


3.9 




AE resolution 


1.7 




Self cross-feed fraction 


0.8 




sPlot bias 


0.6 


Total 




7.4 



dates and estimate an overall tracking reconstruction effi- 
ciency of 2.4%. The uncertainty associated with particle- 
identification (PID) selection criteria is estimated as 1.4% 
by comparing data and Monte Carlo identification effi- 
ciencies for protons from A — > pi: decays and pions from 
Kg — > TTTT decays. The limited signal-Monte Carlo sam- 
ple available to measure the reconstruction efficiency over 
the Dalitz plane results in an additional 0.4% contribu- 
tion. The uncertainty in the number of BB pairs in the 
data sample accounts for a 1.1% contribution, while the 
assumption of a 50% ratio of B^B^ to BB at the r(4S') 
gives an additional 3.2% contribution, computed from the 
difference between 50% and the current measured value 
(48.4 ±0.6)% [10]. 

Uncertainties associated with event-selection require- 
ments on the Fisher discriminant and vertex fit probabil- 
ity are estimated by comparing data and Monte Carlo- 
selection efficiencies for a sample of J/ ipKg candi- 
dates. We use an inclusive sample oi A pir candidates 
to estimate uncertainties associated with the efficiencies 
of the flight-length significance and yl-mass requirements. 

The application of the requirement on the recon- 
structed m(yl7r) invariant mass to suppress B^ — > AcP 
background has two associated systematic effects. The 
first results in an approximate 0.2% increase in the 
branching fraction due to the residual Ac component that 
survives the cut. The second results in an approximate 
0.5% reduction of the branching fraction due to the re- 
duced Dalitz-plot phase space. We correct for these ef- 
fects and take the larger of the two as the uncertainty 
associated with the Ac veto cut. 

We vary parameters that are kept fixed in the like- 
lihood fit by their uncertainties, as measured on the 
signal-Monte Carlo sample fit, and measure the varia- 



tion of the sPlot fitted result. The uncertainties asso- 
ciated with the parameters that enter the definition of 
the signal PDF are conservatively considered as corre- 
lated and are thus added to give a signal-PDF overall 
uncertainty of 3.2%, where the uncertainty in signal- tties 
fixed parameters accounts for a 1.9% contribution and 
that in signal- Ai? fixed parameters for a 1.3% contribu- 
tion. The same procedure is applied to the parameters 
that enter the background PDF, with uncertainties deter- 
mined on luminosity-weighted background-Monte Carlo 
samples, giving an additional 2.2% uncertainty. Finally, 
we combine the two uncertainties in quadrature and ob- 
tain a 3.9% uncertainty associated with the shapes of 
the signal and background models. The comparison of 
B J/ipKg data and Monte Carlo samples reveals that 
the width of the AE Gaussian in the signal PDF can be 
underestimated in the Monte Carlo by up to 5%, which 
translates to an additional 1.7% uncertainty. 

We estimate possible biases associated with the deter- 
mination of yields with the s'Plot technique, using an 
ensemble of Monte Carlo experiments. Signal events, 
generated and reconstructed with a complete detector 
simulation, were mixed with background events, gener- 
ated according to the background PDF. The numbers of 
events were chosen according to the expected yields in the 
data sample under study. We estimate an uncertainty of 
0.6%. 

The main systematic uncertainty in the polarization 
measurement is associated with the limited statistics of 
the Monte Carlo sample used to measure the signal- 
reconstruction efficiency in the (cos^h, E*^) plane, which 
results in Pl {E*{) uncertainties of 0.05, 0.07, and 0.04 
for the three E*^ bins. Variation of parameters fixed 
in the likelihood fit within their uncertainties provides 
additional contributions of 0.004, 0.03, and 0.03 in the 
three bins, respectively. We correct the fit result for 
the small biases we observe in a sample of Monte Carlo 
experiments, where background candidates were gener- 
ated with the helicity distribution observed in tties < 
5.27 GeV/c^ sideband data, and conservatively take these 
shifts as contributions to the systematic uncertainty. 



VI. BRANCHING-FRACTION RESULTS 

We select a total of 6360 candidates in the region 
|A£:| < 100 MeV, tues > 5.2GeV/c^, \m{ATT) - m(Ac)\ > 
20MeV/c2. Table [H reports the fitted values of the 
two-dimensional m^s-AE PDF parameters, while Fig. [2] 
shows projections of the two-dimensional PDF on the 
Toes and AE axes. Figure [3] shows the efficiency- 
corrected signal- s "Plot distribution of candidates as a 
function of m{Ap), demonstrating a near-threshold en- 
hancement similar to that observed in other baryonic B 
decays. Summing the content of the efficiency-corrected 
sVlot bins, we obtain 916 ± 92 signal events, where the 
uncertainty is statistical. Using Eq. [51 we measure the 
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AE(GeV) 



FIG. 2: Upper plot: tties distribution of candidates with 
|A£;| < 27MeV. Lower plot: AE distribution of candi- 
dates with niES > 5.274 GeV/c^. The projections of the 
two-dimensional fit PDF are shown superimposed. 

branching fraction: 

B{B" Apn-) = [3.07 ± 0.31(stat.) ± 0.23(syst.)] x 10"^. 

This measurement, which is compatible with a previous 
measurement by the Belle collaboration 3] , confirms the 
peaking of the baryon-antibaryon mass spectrum near 
threshold, a feature that plays a key role in the explana- 
tion of the larger branching fractions of three-body bary- 
onic B decays compared to two-body decays [H . From the 
maximum-likelihood fit to the branching-fraction asym- 
metry we obtain: 

A = -0.10 ± 0.10(stat.) ± 0.02(syst.), 

which is compatible with zero asymmetry. 

VII. POLARIZATION RESULTS 

Only 3994 candidates populate the E*^ range 
[1.1, 2.4] GeV. Signal candidates are absent in the region 
with E*^ > 2.4 GeV (Fig.[3|) as a kinematical consequence 




2.5 3 3.5 4 4.5 5 

m (GeV/c^) 




1.2 1.4 1.6 1.8 2 2.2 2.4 

E^(GeV) 



FIG. 3: Upper plot: sVlot of the m{Ap) event distribution 
with efficiency corrections applied. Lower plot: sPlot of the 
Ej^ distribution with efficiency corrections applied. Horizon- 
tal bars represent bin ranges. 

of the near-threshold peaking of the baryon-antibaryon 
mass spectrum. 

We plot in Fig. d] the values of the longitudinal polar- 
ization product ctjjP L (-S *^) obtained from the maximum- 
likelihood fit. Table Hn] displays the longitudinal, trans- 
verse, and normal polarization measurements in each of 
the three E*^ bins, assuming = —0.642 ± 0.013 for 
the A decay-asymmetry parameter [lo| . The results are 
consistent with full longitudinal right-polarization of yl's 
from ApTT~ decays at large E*^ (yl's would be op- 

positely polarized). The transverse polarization is not 
expected to be zero because of the presence of strong 
final-state interactions. 



VIII. CONCLUSIONS 

Based on 467 x 10'^ BB pairs collected by the BaBar 
detector at PEP-II, we present a measurement of the 
B'^ — > ApiT~ branching fraction and confirm the peaking 
of the baryon-antibaryon mass spectrum near threshold, 
characteristic of three-body baryonic B decays. In ad- 
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TABLE II: Branching-fraction results. Ns and Nb are the 
numbers of fitted signal and background events, respectively. 
The symbol fi (AE) is the mean for the narrow Gaussian of 
the Ai? signal-PDF component, while ci (A_B) is the slope 
of the linear AE background PDF. /i (m-Bs) is the mean for 
the Gaussian of the rriES signal PDF, and cargus (jties) is 
the coefficient of the exponent in the background rriES Argus 
function [T^]. The uncertainties are statistical. 



Parameter 



Value 



Ns 
Nb 



183.311^:^5 
6176 ± 80 



^^{AE) 
ci (AE) 



-2.65 ± 1.84 MeV 
-3.5±0.4GeV-' 



^i(mEs) 5.2797 ± 0.0003 GeV/c^ 

CARGUS (mEs) -14.6 ± 1.45 




E- (GeV) 



4: The product of A longitudinal polarization and 
Horizontal bars represent bin ranges. 



FIG. 



as a function of E*^ 



dition we measure the A polarization in —> ApTT~ 
decays as a function of E*^. The measurement is com- 
patible with the theoretical prediction of full longitudi- 
nal right-handed polarization for large E*^, which follows 
from the purely left-handed b —^ s transition in the stan- 
dard model [91 . 
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TABLE III; Polarization results. Ns and Nb are the numbers 
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A polarizations in each of the three E*^ bins. 







E\ range ( GeV) 






1.10 - 1.53 


1.53 - 1.80 


1.80 - 2.40 


Ns 


63 ±9 


51 ±9 


55± 11 


Nb 


519 ±23 


643 ± 26 


2663 ± 52 


Pl 


-0.08155:40 ± 0.09 


±0.64lg:^^ ±0.12 


+0-97lo.62 ± 0-08 


Pt 


±0.25lo:58 ±0-09 


+0.56lo:lg ±0.12 


±0.0510.6^ ± 0.08 


Pn 


-QMtlit ± 0.09 


-0.78lo:3e ±0.12 


±0.26lo.53 ± 0.08 
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